o 



Abstract 



Generalised ClifFord groups and simulation 
of associated quantum circuits 

Sean Clarkt, Richard Jozsa^ and Noah Linden§ 

^Department of Computer Science, University of Bristol, 
Merchant Venturers Building, Bristol BS8 1 UB U.K. 
§ Department of Mathematics, University of Bristol, 
University Walk, Bristol BS8 1TW U.K. 

r- 
o 
o . 

Quantum computations that involve only ClifFord operations are classically simulable 
despite the fact that they generate highly entangled states; this is the content of the 
| ^ ■ Gottesman-Knill theorem. Here we isolate the ingredients of the theorem and provide 

generalisations of some of them with the aim of idcntifying new classes of simulable 
quantum computations. In the usual construction, Clifford operations arise as projective 
normalisers of the first and second tensor powers of the Pauli group. We consider 
replacing the Pauli group by an arbitrary finite subgroup G of U(d). In particular 
we seek G such that G ® G has an entangling normaliser. Via a generalisation of the 
Gottesman-Knill theorem the resulting normalisers lead to classes of quantum circuits 
that can be classically emciently simulated. For the qubit case d = 2 we exhaustively 
treat all finite subgroups of U(2) and find that the only ones (up to unitary equivalence 
\ and trivial phase extensions) with entangling normalisers are the groups generated by 

C-». ■ X and the n th root of Z for n £ N. 

o 

Ph! 1 Introduction 

The identification of classes of quantum computations that can be classically emciently 
simulated is a bàsic tool for studying the relationship between classical and quantum com- 
putational power. One of the earliest results in this context is the Gottesman-Knill (GK) 
theorem [H [2] • It initially arose in the study of the stabiliser formalism for quantum error 
correcting codes and has a rich variety of mathematical ingredients. In this paper we isolate 
these ingredients and develop generalisations of some of them with an aim of identifying 
further new classes of simulable quantum computations. (Other directions of generalisation 
of the GK theorem were considered in [5]). 

The study of classical simulation is closely related to the invention of novel mathemati- 
cal formalisms for the representation of quantum states and computational steps, and the 
conseqüent calculation of measurement probabilities. Indeed relative to any such mathe- 
matical formalism there will be a class of states with "small" descriptions whose updates 
can be emciently calculated within that representation. The most commonly used formal- 
ism describes states in terms of their amplitudes in the computational basis and then the 
property of entanglement leads generically to exponentially large descriptions, inhibiting 
emeient classical simulation. Consequently[3] we may claim that entanglement is an im- 
portant resource for quantum computational power in the sense that if it is absent then 
the quantum computation may be classically emciently simulated. However the situation 
becomes less clear when we realise that in alternative formalisms the class of states with 
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suitably small descriptions can include rich varieties of entangled states leading to the efi- 
cient simulation of classes of computations that involve much entanglement along the way. 
In the development of the theory of quantum error correction, highly entangled quantum 
code states appeared increasingly ungainly in the Standard amplitude description. The 
stabiliser formalism was introduced pQ to provide a compact and manageable description 
leading to the notion of general stabiliser states, having small descriptions in this formal- 
ism. Subsequently one could ask for quantum operations that preserve this class of states 
and hence arrive at the GK theorem, loosely speaking, that associated classes of quantum 
circuits (although developing entanglement) can be efficiently simulated. 

To motivate our proposed generalisations we first outline the ingredients of the GK 
theorem. Let X and Z denote the Standard qubit Pauli operations. The Pauli group is 
defined by T = {X, Z, il) (where the pointed brackets denote the group generated by the 
enclosed elements). The Pauli group on n qubits is defined as the n-fold tensor power 
T n = T® n which is a finite group of size \T n \ = 0(4 n ). If \tp) is any n qubit state define its 
stabiliser as 



Clearly Stab(^>) is always a subgroup of T n (albeit the trivial subgroup for many \ip) ) s). 
IV') is a stabiliser state if it is uniquely characterised by its stabiliser i.e. it is the only state 
left invariant by all g G Stab(^>). Any subgroup of a finite group G may be described by 
0(log |G|) generating elementsp], providing our small description of stabiliser states (which 
can generally be highly entangled): Stab(^) = (gi, ■ ■ ■ ,g r ) with r = 0(n). 

Next note that for any U, g \tp) = iff (UgU^)U = U so corresponding to the 
application of a gate U to \ip), Stab(^) is updated by conjugation of the generators by U. 
In order to guarantee that U \ is again a stabiliser state we require that UgUÏ G T n for 
all generators g. To guarantee this condition independently of the choice of we impose 
it for every g G T n introducing the so-called Clifford operations: a Clifford operation on n 
qubits is a unitary n-qubit operator C with the property that CgC^ G V n for all g G V n . 
For each n we obtain the Clifford group C n defined by 



i.e. C n is the (group-theoretic) normaliser of the group T n (within the unitary group). Let 
H denote the Hadamard operation, let P denote the 7r/4-phase gate 



and let CZ denote the 2-qubit controlled-Z gate. Then we have the following full explicit 
characterisation [U [T] of Clifford operations. 

Lemma 1. C\ = (H, P) and C2 = {C\ ® Ci, CZ). For n > 3 an n-qubit gate U is in C n iff 
it can be expressed as a circuit of gates from C\ and C2 ■ 

In terms of these structures we can give a precise statement of the GK theorem: 

Theorem 1. (Gottesman-Knill theorem). Consider any polynomial-tíme quantum compu- 
tatíon of the following sort. The starting state is a stabiliser state and each computa- 
tional step is one of the following: 



Stab(V) = {g G T n : g Wl = |V>}}- 



C n = {C G U{2 n ) : CTn& = T n } 
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(a) a measurement on a qubit in the Z basis, or 

(b) applícation of a gate from C\ or C2 (not depending on measurement outcomes in (a)), 
or 

(c) applicatíon of a gate from C\ or C2 chosen adaptively depending on a previous measure- 
ment outcome from (a). 

Finally the output is the result of a measurement of the first qubit in the Z basis. 
Then the computation may be classically efficiently simulated. 

The Standard proof of this result (see e.g. [2] or [5] for a more recent improved algo- 
rithm) proceeds by updating the stabiliser description of the state through the course of 
the computation. The update procedure for (b) (and (c) once the measurement result is 
given) is via the normalising property of the Clifford group in relation to the Pauli group 
V n containing the generators. This purely group-theoretic relationship in itself, may be 
entertained for any group G replacing the Pauli group. On the other hand the stabiliser 
update rules for (a) (as elaborated for example in [2j page 463) depend on further features 
specific to the Pauli group, such as the fact that in this group every two elements either 
commute of anti-commute. 

The starting point for our generalisations is an alternative simpler proof in the absence 
of the adaptively chosen steps in (c) : instead of forwardly propagating the state description 
we will backwardly propagate the final measurement allowing us in particular even to free 
the simulation from requiring stabiliser states. Thus let C now be any circuit of Clifford 
operations on starting state |V>o) which is now not required to be a stabiliser state. If the 
final measurement on the first qubit has outputs 0,1 with probabilities po,pi then po — p\ 
is given by the expectation value of Z\ = Z (g) / (g) . . . (g> / in the final state C 

Po-PiH^lc-tZiChfo). (i) 

This computation suffices to simulate the output (as we also have po+pi = 1). Now Z\ £ V n 
so C^Z\C has the product form Pj x <g) . . . ® P, ln for Pauli operators Pi k . Hence if |Vo) 18 an y 
product state |V>o) = \o>i) ■ ■ ■ \&n) then we get 

n 

Po -Pl = JJ (o-k\Pi k Wk) (2) 
k=l 

which can clearly be calculated classically in linear time 0(n). Similarly the commuting of 
the successive one and two qubit Clifford gates through Z\ can also be done in time linear in 
the size of the circuit giving a linear time classical simulation of the quantum computation's 
output. 

This approach to the simulation of Clifford circuits may also be extended to allow for 
measurement steps (of type (a) above) so long as subsequent gates are not chosen adaptively 
(as they are in (c) above, with stabiliser starting states). To achieve this we replace each 
measurement step by the following: for each measurement on a qubit i adjoin an extra 
initial qubit in state |0) and replace the measurement step by a (Clifford) CNOT operation 
with control and target being the i th and new qubits respectively. The newly introduced 
qubit is not used in any other way by the computation so its presence serves to decohere 
the z th qubit into the post-measurement mixture i.e. each measurement step of the form 
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(a) is replaced by a CNOT step of the form (b) and the final output is unchanged. We may 
ask if a further such trick could allow efncient simulation of the output of the process with 
adaptively chosen gates (as in (c)) in addition to just measurements (a) themselves, for the 
scenario of Clifford circuits on arbitrary product starting states. Such further generalisation 
is not likely to be possible for the following reason: if we allow arbitrary product state inputs 
and adaptive Clifford gate choices then we could (as shown in [6j) implement the 7r/8-phase 
gate 

S = ( e^ 4 ) ' 

To see how this is achieved let \a) = 4j(|0) + e í7r//4 |1). Then for any qubit \if>) ap- 
ply CNOT to \ip) \a) and measure the second qubit. If the outcome is then the post- 
measurement state is S \ip) |0). If the outcome is 1 then the post-measurement state is 
£t |1) e l?r//4 so applying P to the first qubit gives S \vp) |1) up to an overall phase. Thus 
we implement S in either case by responding adaptively to the measurement outcome. A 
supply of \a) states can be provided as an extension of the input product state. Now it is 
known that S together with C2 is a universal set of gates for quantum computation so we 
would get an efhcient simulation of all poly-time quantum computation, which is generally 
believed not to be possible. 

Our discussion above generalises the GK theorem by allowing arbitrary product state 
inputs (and noting that also arbitrary entangled stabiliser state inputs can be generated by 
a prefixed Clifford circuit on the product state |0) . . . |0)) but on the other hand restricts 
the original form by not allowing adaptive choices of gates. Its virtue is that it relies only 
on the group-theoretic normaliser relationship between Clifford and Pauli groups and may 
thus be immediately generalised to having arbitrary unitary matrix groups G replacing the 
Pauli group as the starting point. We require no associated subgroup structure to support 
a stabiliser state formalism nor any consideration of stabiliser states themselves. 

Let G be any finite matrix subgroup of U{d). Introduce the (linear) normalisers of G 
and G®G: 

N{G) = {U £ U(d) : UGW = G} 
J\Í(G ®G) = {U e U(d 2 ) : U(G íg> G)rf = G®G}. 

We will be interested in using normaliser operations as circuit gates and we can therefore 
allow extra overall phases to be generated in the above relations. Thus we introduce the 
notion of projective normaliser. 

VM(G) = {U G U(d) :VgeG, UgW = cg' for some g' G G and c G §1} 
VAÍ(G ®G) = {U G U{d 2 ) :VgeG®G, UgU^ = cg' for some g' G G <g> G and c G §1} . 

where §1 = {é e : < 6 < 2vr}. 

Remark. The significance of projective normalisers is illustrated by the following example. 
If G is the Pauli group (X, Z, il) then N(G) = VM{G) and it contains the phase gate P. 
But if G is the group (X, Z) comprising matrices with only real entries then M(G) does 
not contain the phase gate P but VM{G) does capture this gate remedying the absence of 
complex elements in the centre of this smaller real number version of the Pauli group. □ 

Mimicking our previous discussion we will especially seek examples of groups G such 
that VM{G (g> G) contains an entangling gate (such as CNOT in the case of G being the 
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Pauli group). Otherwise all normaliser circuits will preserve product states and be com- 
putationally uninteresting. Furthermore in view of eq. ([T]) it is desirable that G contains 
a Hermitian element (such as Z) which can be associated with a measurement. Then we 
will be able to efiiciently calculate its expectation value in the final state of any normaliser 
circuit (with product state input) forming the basis of our classical simulation procedure. 
Even if G does not contain a Hermitian element we may use the Hermitian matrix A + 
for any A £ G and similarly apply the arguments following eq. (PQ) to simulate an associated 
measurement expectation value. In this vein we also note that if G acts irreducibly on C d 
(e.g. as is the case for the usual Pauli qubit group) then any d x d matrix may be expressed 
as a linear combination of the matrices of G (c.f. [18] p. 48) so we may emciently com- 
pute the expectation value for any von Neumann measurement on a single qudit or more 
generally on O(logn) qudits. 

1.1 Teleportation groups 

In this paper we will consider only subgroups G of U(d) that act irreducibly on C d . In ad- 
dition to facilitating the mathematical analysis at various stages such groups have an extra 
significance as prospective generalised substitutes for the Pauli group as follows. Recali that 
another fundamental appearance of the Pauli group is in quantum teleportation, providing 
the set of Bob's "correction operators". In measurement based quantum computation[7], 
which can be viewed from the perspective of teleportation [8] the associated Clifford oper- 
ations have a special role of being parallelisable to depth 1 in this formalism. Thus we may 
ask what other sets of operators may appear as Bob's correction operators in generalised 
teleportation schemes and then ask for their normalisers. This will again lead to classes of 
computations that are parallelisable in the corresponding generalised measurement based 
computational model. In this regard, irreducible subgroups of U(d) play an important role. 

Let us define a generalised teleportation scheme as follows. Alice and Bob share the 
maximally entangled 2-qudit state |0) = — ^ \ï) \i) and Alice also has a 1-qudit state 
|a). Let A4 = {A±, . . . , A r } be any 2-qudit generalised measurement (POVM). Suppose 
Alice applies Ai to the first two qudits of \a) 1 |çí>) 23 obtaining measurement outcome Ai. 
Let pi be Bob's post-measurement state. (We may without loss of generality take the full 
post-measurement state to be \J~A~i \a) \cp) renormalised, and pi is obtained by tracing out 
the first two qudits). This comprises a generalised teleportation scheme if there exists a 
family of unitary operators Ui parameterised by the measurement outcomes, such that for 
all |a) and all i, pi is the pure state Ui \a) i.e. U\ functions as Bob's correction operator 
for measurement outcome i. In the case that {U\, . . . , U r } also forms a group we have the 
following. 

Lemma 2. Let G = {U\, . . . , U r } be any finite subgroup of U{d) that acts irreducibly on 
C d . Then there exists a generalised teleportation scheme with G comprising Bob's correction 
operators. 

Proof: Define \a.- t ) = U\ (g> / |(/>) for i = 1, . . . , r and introduce the positive rank 1 operators 
Ai = ^ \ ai) (ai\ . Then using Schur's lemma (by virtue of the irreducible action of G) we can 
see that ^ i Ai = I d 2 so {Ai, . . . , A r } is a (rank 1) POVM. Furthermore a straightforward 
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calculation gives 

V^ï|«)l |0>23 = -TjTf I a i)l2^l«>3 
V l G l 

and thus Bob's post-measurement state is U \a) as required. Also each measurement out- 
come occurs with equal probability 1/|G|. □ 

In view of this result we introduce the term teleportation group to refer to any finite 
subgroup of U(d) that acts irreducibly on C d . 



1.2 Outline of the paper 

Returning to our primary motivation of classical simulation we would ideally wish to find 
all teleportation groups G in U(d), compute the projective normalisers of G and G <S> G 
seeking especially the cases of G such that VAÍ(G <g> G) contains an entangling gate. We 
refer to such groups as entangling teleportation groups. 

In the case of the Pauli group the projective normalisers are known explicitly ana- 
lytically. However the derivation is lengthy and rests on many properties special to the 
Pauli operators. We are not able to similarly explicitly analytically characterise projective 
normalisers for general teleportation groups (even for d = 2) and we resort to exhaustive 
methods using various computer àlgebra packages. In the qubit case d = 2 we will be able 
to treat exhaustively all possible teleportation groups. In section [2] we will describe our 
algorithm for computing normalisers and projective normalisers of G and G (g) G for any 
given teleportation group. In these methods it will be important to cut down wherever 
conveniently possible, the range of various cases that needs to be considered to allow the 
computer àlgebra to terminate in a reasonable time. In this respect it is important to note 
that the centre Z{G) of any teleportation group G (which by Schur's lemma comprises only 
phase múltiples of the identity) plays no role in extending or limiting the existence of pro- 
jective normalisers. Hence in section [3] we will describe how a search over all teleportation 
groups in U (d) for entangling ones, can be reduced to the study of normalisers of projec- 
tively inequivalent projective representations of the so-called base groups in U(d) which are 
defined to be the central quotients of teleportation groups. 

Next, in section [4] we will apply our methods to identify all possible entangling telepor- 
tation groups in the qubit case of subgroups of U(2). We prové the following result. 

Theorem 2. The only finite subgroups G of U{2) (up to unitary equivalence and trivial 
phase extensions) such that G®G has an entangling projective normaliser, are (X, Z 1 /* 1 ) 
for n £ N (the usual Pauli group being a central extension of the case n = 1). 

Finally in section [5] we will make some concluding remarks and identify some avenues 
for further developments. 



2 Algorithm for determining normalisers and projective nor- 
malisers 

In this section we describe a procedure for computing the linear and projective normaliser 
elements of a teleportation group G = {Uj} C U(d) 
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Let Gen(G) C G be a set of generators for G. For each U £ Gen(G) let U' £ G denote 
the image of Í7 under conjugation with some N £ N(G): 

NUN ] = U'. (3) 

Let us rewrite the normaliser matrix TV as a (cí 2 x l)-column vector n where 

^dj+fe = N jtk . (4) 

Then eq. ([3]) becomes 

(J (8) C/ T - U' I)n = (5) 
where I is the (d x d) identity matrix. 

By specifying the vàlues of one such pair U and U 1 and treating the entries in the vector 
n as unknowns we can obtain from eq. ([5]) d 2 simultaneous equations in d 2 unknowns. If we 
assign members of G as the images of all the elements of Gen(G) then we can solve these 
equations simultaneously by finding the null space of (I®U r -U'®I) for each U £ Gen{G) 
and its chosen image U'. If a non=trivial solution exists for n simultaneously for all U then 
this gives us a solution for TV £ J\Í{Q). 

This provides us with the basis for an algorithmic approach to computing the elements of 
the normaliser. We enumerate all the possible choices of images of the generators of G and 
solve the simultaneous equations, discarding the trivial solutions. In order to improve the 
performance of this approach we observe that all the mappings on G induced by conjugation 
with a normaliser element ./V are constrained by the fact that each image must have the 
same order as the generator and any pairwise choice of images must preserve the group 
relations of the corresponding generators. 

Thus we get the following algorithm for computing the normaliser elements of a tele- 
portation group. 

Procedure 1 - Compute M(G) 

1. For each U in G compute Order(U), the elements of G with the same order as U, 
Comm(U), the elements of G that commute with U and NComm(U), the elements 
of G that do not commute with U. 

2. Take a minimal set of generators Gen{G) = {Í7i, . . . U r } of G and find the set of all 
pairs that commute. 

3. Calculate all possible images of Gen(G) by considering all the choices given in steps 
4,5 and 6. 

4. For the images of the first generator, U\, iterate through the set Order(Ui). 

5. The set of possible images of each subsequent generator Uj is formed by starting with 
the set Order(Uj) and then repeatedly intersecting with Comm{Uk) if Uj and Uf. 
commute and with NComm(Uk) otherwise for each k < j. 

6. For each choice of possible images {U[, . . . U^,} of {U\, . . . U r } compute the combined 
null space of (J ® Uj - U'j (8) J) for j £ {1, ... r}. Any non trivial solution corresponds 
to a normaliser element of G. 
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We can apply the same procedure to compute N(G (8> G). In addition we may also 
test to see if a normaliser gate is entangling using the following result. A 2-qudit unitary 
operator V G U(d 2 ) is said to be entangling if for all A,B<E U(d) it is true that V 7^ A® B 
and V / SWAP(A®B) (where the SW AP operation is defined by SWAP \i) \j) = \j) 
Then we have the following characterisation [9]: V is not entangling if and only if one of 
the two following conditions holds for every i,j, k, l, i, j ,k,l G {0, . . . ,d — 1}: 

By checking these simple algebraic conditions we can readily identify if a given operation is 
entangling or not. 

2.1 Algorithm for projective normaliser elements 

To develop an algorithm for determining projective normalisers of a group G we first show 
that any such element can be found as a linear normaliser of a group G' generated by adding 
suitable additional central elements to G. 

Let N be any (fixed, chosen) projective normaliser element for G. Then for all Í7 G G 



Since any U G G has C/' G ' = I, c must be a \G\ th root of unity. Thus if G 1 is the group 
obtained by including all such roots into G we see that any operator is a projective nor- 
maliser of G iff it is a linear normaliser of G' . In practice (especially when treating larger 
groups such as G (g> G) this extension of G to G' becomes too large to be manageable for 
subsequent application of exhaustive enumerations in procedure 1. Thus we develop more 
refined restrictions on c to further limit its possible vàlues. 

Note first that there is ambiguity in the choice of c and V in eq. (|6|) due to central 
phases that may already exist in G. This is remedied using the following lemma. 

Lemma 3. If G is any teleportation group then the centres of G and G®G are both cyclic, 
comprísing phase múltiples of I. 

Proof. Since G acts irreducibly Schur's lemma guarantees that any central element is 
a múltiple of /. Thus Z(G) is a finite subgroup of §1 and hence is cyclic (necessarily 
generated by its element e %e with least positive 0.) For G (g> G let g± (8> 52 be any central 
element. Then it commutes with g ® g^ 1 for all g G G so g\ G Z(G). Similarly 52 G Z(G) 
so Z(G ® G) = Z(G) 8) Z(G) is again a finite subgroup of Si, hence cyclic. □ 

Now let oj s I with uj s = e 2m / s be the minimal phase element of Z(G). Then in eq. ([6]) 
we choose c = e ld such that < 6 < 2ir/ s which fixes c and V uniquely. Furthermore with 
this choice, the unique correspondence between U and V means that we can view c as a 
function of V (rather than U): 



NUNÏ = cV with c G §1 and V G G 



(6) 



NU = f{V)V with < arg{f{V)) < 2vr/s. 



(7) 
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(The function / will also depend on the choice of N but we omit explícit inclusion of this 
parameter for notational clarity.) We will also refer to the association of phase vàlues /(V) 
to V G G as a phase function for G. Introduce 

r = {f(V)V -.V eG} = NGNl 

Thus r is a unitary matrix group isomorphic to G and Z(T) = Z(G). 

Lemma 4. Let {U\, . . . , U r } be any set of generators for G. 

(a) Then {uj s , f(Ui)Ui, . . . , f(U r )U r } generatesT. 

(b) If Ui has order ni then f(Ui) has the form LOg n . for some < k < sn«. 

Proof. (a) For any V G G we have V = U h U Í2 ...U im = \\ k U ik . Also from eq. (0) 
/ : G — ► Si has the multiplicative property: /(V1V2) = f(Vi)f(V2)z for some z G Z(G). 
Hence f(V)V = z\[ k f{U ik )U ik for some z G Z(G). Thus cu s I together with f{Ui)Ui for 
i = 1, . . . , r generates T. 

(b) We have f(Ui) n 'U^ = f(Ui) n *I which is thus in Z(T) = Z(G). Hence /(t/i)™ 1 = ^ so 
/(i7i) is a power of w„ íS . □ 

For any projective normaliser of G, lemma 0] provides restrictions on the vàlues that 
the phases /(V) in eq. ([7]) can possibly take. Define <J?(G) to be the set of all choices of 
f(U\), . . . , f(U r ) satisfying the conditions (a) and (b) of the lemma. Extending G by new 
central elements $(G) will then give a group G' whose linear normalisers are precisely the 
projective normalisers of G. In many practical examples |Z(G)| = s and the generator 
orders rtj are small compared to |G| so the resulting extension to G 1 can be far smaller than 
that obtained by simply adding all |G| th roots of unity to G. Correspondingly we introduce 
the following computational procedures. 

Procedure 2 - Compute <&{G) 

1. Take a set of generators {Ui, . . . ,U r } of G. Let the orders of the generators be 
{m, . . . ,n r }. 

2. Take a generating element z of the centre of G. 

3. For each possible combination of jt G {0, . . . , n^s — 1} for t = 1, . . . , r perform steps 4 
and 5. 

4. Let f(U t ) = uÍ\ s . 

5. If (z, f{U t )U t : í G {1, . . . ,r}) = G then add each f{U t ) to the set $(G). 

6. Output $(G). 

We note that, for small groups, step 5 can be performed relatively quickly using a 
computational package such as GAP[15j. 

This provides us with an algorithm that produces a group G' such that the projec- 
tive normalisers of G are the linear normalisers of G 1 and hence a means to compute the 
projective normaliser elements of a teleportation group and its tensor square. 

Procedure 3 - Compute VJ\Í(G) 
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1. Compute ^(G) using procedure 2. 

2. Compute G' = (<j>I, U : <fi £ $(G), U G G). 

3. Compute and output the linear normaliser elements of G' using procedure 1. 

Procedure 3 can be used to compute if teleportation groups are projectively entangling 
in the following manner. We compute the two-qudit matrix group G <S> G from G and use 
procedure 3 applied to the group G ® G to find VJ\f{G (8> G) (recalling that the centre of 
G®G is cyclic). Each projective normaliser element found can then be tested to see if it is 
entangling. 

3 Base groups and their projective representations 

Our principal aim is to apply the procedures of the preceding section to systematically 
study classes of teleportation groups. In the next section we will exhaustively treat the 
qubit case of all teleportation subgroups of Z7(2). There are infinitely many teleportation 
groups in U (d) but (at least for small d) they are known to fall into regular families. Instead 
of directly enumerating these we will adopt a different approach with a view to reducing 
the amount of computer àlgebra required. It is clear from the definition of VN{G) that the 
group's centre Z{G) plays no role in restricting or enabling new projective normalisers. 

For any teleportation group G we introduce the central quotient B = G/Z(G). These 
central quotients are called base groups of U{d). (Our term "base group" is synonymous 
with "finite collineation group" in pH]). Let T : B —* G be any chosen transversal of 
Z(G) in G i.e. a choice of element in each coset of Z{G) in G. We also require that 
I is chosen from Z(G) itself. By slight abuse of notation we will also use T to denote 
the set of matrices {T(b) : b G B}. Thus T defines a projective representation of B 
i.e. for all 61,62 G B T(b\)T(b2) = cT(6i62) for some c G Si. Conversely any projective 
representation p of B arises as a transversal of a representation of some G with B = G/Z{G). 
(Indeed G may be generated by the matrices of p together with extra central phases c from 
p(6i)p(6 2 ) = cp(bib 2 )). 

Lemma 5. Let T be any projective representation of B = G/Z{G). Then N is a projective 
normaliser of G (resp. G<S>G) iff N is a projective normaliser of the set T (resp. T®T = 
{A®B:A,BeT}). 

Proof: immediate from the fact that every U £ G has the form U = cV for some V G T 
and c G §1 (and similarly for G ®G and T ® T). □ 

Let Gi and G2 be teleportation groups with isomorphic central quotients B. Choose 
transversals T\ and T2 giving projective representations of the central quotients. We say 
that Gi and G2 are projectively equivalent if the projective representations Ti,T 2 (for some 
hence any choice of transversals) are projectively equivalent as projective representations 
of B i.e. there is A G U(d) and c(6) G §1 such that ï\(ò) = c{b)AT 2 {b)A^ for all 6 G B. 

Recali that G is called entangling if VN{G (8> G) contains an entangling operation. 

Lemma 6. Let G\, G 2 be projectively equivalent teleportation groups. Then G\ is entangling 
iff G 2 is entangling. 
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Proof: Let Ti,T2,A be as above. Suppose ./V is an entangling projective normaliser for 
Gi. Introduce M = {A® A)N{A^ <g> A^). Since Ti is a transversal of the centre of Gi every 
member U of G{ has the form cV for V G Tj and c E §1. Thus a straightforward calculation 
shows that M is a projective normaliser for G2 ® G2. Also M is locally equivalent to iV so 
it is also entangling. □ 

In view of the above lemmas, to find all entangling teleportation groups G it suffices 
to look at a complete set of all projectively inequivalent projective representations T of all 
base groups B and determine if the projective normaliser of T <8> T contains an entangling 
operation or not. G up to unitary equivalence is then generated by the matrices of T 
and further central phases cl. Actually even a complete list of projectively inequivalent 
projective representations may involve redundàncies as the normaliser structure we seek is 
a property of a set of matrices irrespective of how the set represents a group. Then note 
that it is possible for two (projective) representations of a group B to be inequivalent yet 
comprise the same overall set of matrices (up to overall phases) which are then associated 
with the elements of B in different ways (c.f eq. (|24p later for a non-projective example). 

To access the full list of base groups for d = 2 (and also for d = 3, 4) we note that 
the full list of subgroups of SU(2), SU(3), SU(4) are known and provided in [ÏD] (with the 
latter two cases considered in more detail in [TT] and [12] respectively.) 
Remark. These lists are complete up to abstract group isomorphism but not complete up 
to unitary equivalence. In this regard it is important to note that the normaliser group 
M(G) is not a property of an abstract group but of a given representation i.e. two unitarily 
inequivalent representations of the same group will generally have different normalisers. □ 

To pass from base groups of SU(d) to those of U(d) we have the following. 

Lemma 7. The sets of central quotients of finite subgroups of SU(d) and U{d) are the 
same. 

Proof: If G Ç SU (d) then G Ç U(d) so its central quotient is in both sets. Conversely if 
G Ç U(d) then the (finite) group G' generated by det(g)^ 1 ^ d g for all g G G (and any choice 
of d th root) has G' Ç SU(d) and the same central quotient as G i.e. G/Z{G) also appears 
in both lists. □ 

Hence the full list of base groups of U{2) (as abstract groups) is obtained from the 
central quotients of the lists in [TÜ] . 

In order to find all the projectively inequivalent irreducible projective representations of 
a finite group we may use the concept of a covering group p~H [Ï3] and it sufnces to calculate 
the inequivalent irreducible linear representations of this group. 

Definition 1. (page 361 A covering grou^ G* of a finite group G is a finite group 

which is an extension of G with kernel contained in the centre of G* such that every projec- 
tive representation of G is equivalent to one which can be lifted to a linear representation 
ofG*. 

In this construction the linear representations of G* arise from the projective represen- 
tations of G by inclusion of further central elements (phase múltiples of the identity). 

1 Both 'covering group' [TJ] and 'representation-group' [13 are used in the literature as translations of 
Schur's 'Darstellungsgruppe'. We prefer 'covering group' so as to avoid the confusion of constructing a group 
representation of the representation-group. 
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The existence of such groups is then established by the following theorem. 

Theorem 3. Every finite group G of order n has at least one covering group of order nm 
where m is the size of the Schur multiplier (2nd cohomology group) of G. 

Proof: Originally due to Schur. See Karpilovksy [T3] , 

Corollary 1. In each dimension d>2 there is a finite number of projectively inequivalent 
irreducible unitary projective representations of each finite group G. Each of these can be 
lifted to an irreducible unitary linear representation of a covering group of G. 

Proof: This follows directly from theorem [3] specialised to irreducible unitary represen- 
tations of a particular dimension. 

Now we are in the position, for a given dimension d, to find all the inequivalent telepor- 
tation groups using the following procedure. 

Procedure 4 - Find projectively inequivalent teleportation groups in U(d) 

1. Let S be a set of finite subgroups of SU(d) up to isomorphism. 

2. Calculate all base groups of U(d) as the set of central quotients of the elements of S. 

3. Let Bc denote a set of covering groups for the base groups. 

4. A complete set of projectively inequivalent teleportation groups of U(d) is generated 
by the set of inequivalent irreducible unitary linear representations of the elements of 
B c - 

We now provide some additional practical information for the steps in procedure 4. 

1. The lists of finite subgroups of SU(d) for d = 2,3,4 are provided in the literature. 

2. In practice central quotients may be calculated using a computational package such as 
GAP [15] . The equivalence of the base groups of SU(d) and U(d) is shown in lemma 

m 

3. The existence of the covering groups in step 3 comes from theorem [3] and in practice 
we use GAP to calculate them. 

4. The set of inequivalent irreducible unitary representations of the covering groups can 
be produced in GAP. 

5. For each resulting representation G we apply procedures 1,2,3 to determine whether 
G (g> G has an entangling projective normaliser or not. 

3.1 Teleportation groups represented in GL(d,C) 

In applying the above procedures we perform computations using the GAP computational 
system [15] and in particular the REPSN [Ï6] representation theory package to find the 
irreducible representations of finite groups. For some groups REPSN does not produce a 
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unitary representation but a general linear representation in GL{d, C). In this section we 
relate the previous results concerning the normaliser and projective normaliser of faithful 
irreducible unitary representations to tl·ieir general linear counterparts and show that the 
algorithms we provide may be performed with faithful irreducible representations of a flnite 
group in GL(d, C) and the results applied to the faithful irreducible unitary representations. 

Two general linear representation p L , p\: G — > GL(d, C) of a finite group G are said to 
be equivalent if there exists M E GL(d, C) such that for all g E G 

p 1 L (g) = Mp 2 L (g)M- 1 . (8) 

Theorem 4. For any finite group G and faithful irreducible representation pi : G — > 
GL(d,C) of G there exists a faithful irreducible unitary representation pu : G —>U(d) and 
a matrix E E GL(d, C) such that for all g E G 

Pu {g)=Ep L (g)E-\ (9) 

Proof: see [UJ page 74. □ 

Furthermore if pi in the theorem ranges over a full list of inequivalent linear represen- 
tations then pu ranges over a full list of unitarily inequivalent representations. For general 
linear representations we can define linear and projective normalisers as follows. 

^GL(d,C)(PL) = {Ne GL(d, C) : Npi(g)N~ 1 E p L (G) for all g E G} . (10) 

VM GL{dtC) ( PL ) = {N E GL(d,C) : Vg E G 3g' E G,c E C : Np^N' 1 = cp L (g')} (11) 
Theorem H] immediately gives: 

Corollary 2. Given the two representations pi and pu of theorem^ then every normaliser 
element N E J\fu{d)(pu) defines a normaliser element M = E~ l NE E Ngl^qÍPl)- 

Note that in eq. (jl lfí c must actually be in §i since G is finite. Thus our previous pro- 
cedures 1,2,3 can be used unchanged to compute the normalisers and projective normalisers 
of pl and pl ® p^. If the latter fails to be entangling we can conclude by the corollary 
that pu is also not entangling, without having to carry out the translation from pl to pu 
explicitly. 

4 Application to teleportation groups in U(2) 

We now apply our procedures to consider all teleportation groups in U(2) and show that the 
entangling ones are exactly those which have a central quotient isomorphic to a dihedral 
group of order 4m for some integer m. Alternatively these teleportation groups may be 
described as unitary equivalents of groups obtained by adjoining additional central elements 
to the matrix groups (X, Z 1 /' 11 ) for n E N. 

We use a labelling system for small finite groups that is used in the GAP computational 
package. This system assigns two numbers to a group. The first is the order of the group 
and the second is a unique index for each group of a particular order. As examples we have 
[4, 1] for the cyclic group of order 4, [4, 2] for the Klein four group, [12, 3] for the alternating 
group on 4 elements etc. 
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4.1 Base groups of £7(2) 



We use the list of base groups of SU{2) which are given in Blichfeldt [10] as the finite 
collineation groups and which we are also the base groups of U{2). These consist of two 
infinite families of groups and three 'special' groups. The infinite families are the cyclic 
groups, which have no irreducible representations in U(2), and the dihedral groups Di n of 
order 2n. The three special groups are the tetrahedral group A4 = [12,3], the octahedral 
(or cube) group S4 = [24,12] and the dodecahedral (or icosahedral) group A§ = [60,5]. 
First we deal with the three special groups and fmally the dihedral case. 



4.2 The tetrahedral group as base group 

A covering group of the tetrahedral group A4 = [12,3] is [24,3] [13] . We use the compu- 
tational package GAP to compute three inequivalent representations in U{2) which is the 
maximum number of inequivalent representations. These are the matrix groups M\, M2 
and M3 where 

Mi = 
M 2 = 
M 3 = 

Using an implementation of the algorithm to find normalisers given in section [2] it can 
be seen that M\ , Mi and M3 are not entangling. To see that these three representations are 
not projectively entangling the ranges of the possible phase functions of the matrix groups 
Mi, M2 and M3 are calculated using the algorithm described in section ETTl From this we 
find that the possible phase functions take vàlues in |cc^ : j = 0, 1, 2^. When we add the 
central elements corresponding to these phases to the matrix groups Mi,M2 and M3 we 
get matrix groups isomorphic to Zq ® [24, 3] = [72, 25]. This group has only one irreducible 
representation in U(2) up to equivalence. This can be represented as 

(US f) > (7 3)) w 

and we compute that it is not entangling. Hence no teleportation group in U(2) with central 
quotient isomorphic to the tetrahedral group is entangling. 
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4.3 The octahedral group as base group 

The octahedral group S4 = [24, 12] has [48, 29] as a covering group. This has one rep- 
resentation in U(2) up to equivalence which is not entangling. The generators we used 
are 



V2 v -V'V2W A)l 



(16) 
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All possible phase functions take vàlues in {1, i} so to test to see if this group is projectively 
entangling we must test if any linear representation of [48, 29] = [96, 192] is entangling. 
Up to equivalence there is only one faithful representation in U(2) and that has generators 

TiO -O'Titi 3))- (17) 

We have computed that this matrix group is not entangling. We conclude that no telepor- 
tation group with central quotient isomorphic to the octahedral group is entangling. 



4.4 The dodecahedral group as base group 

The dodecahedral group A§ = [60,5] has unique covering group [120,5]. There are two 
faithful irreducible representations of this group up to equivalence in U(2). The computa- 
tional system gap gives us representations in GL(2, C) which by the results of section 13.11 
will suffice for our calculations. The first representation has the following two generators 



, .2 1 , ,4 , .8 , ,11 , ,14 9, , 9, ,4 , ,7 , ,13 

W15 - " í5 + u l5 - u l5 - u 15 - u 15 -2u í5 - 2uj 15 - u 15 - u 15 

, ill ,,14 , , 4 , ,7 , ,,11 ,,13 i,, 14 

-u í5 - u 15 -ÜJ Í5 - u 15 - u l5 + u 15 - u í5 + u l5 

-U 15 + uf 5 - uf 5 + uf 5 + u\\ + u\l UX5 + uf 5 

"h + < + Ml + 2 "ll "15 + tvf 5 + u\ 5 - u\l + "II - "ll 

The second representation has the following two generators 

1 /-2W15 - 2uf 5 - U 7 15 - u\l ~ "II - "ll M 5 + "15 + M 5 + "ll + ^15 + "ll 

2 V ~"ï 5 + "ll " "ll + "lí "l5 " "15 + ^15 " "ll 



(18) 

(19) 
(20) 



1 ("15 + "h + "l5 + "L + ^15 + 2^15 + W ÏI + 2 "\í "15 ~ "h + "íï + "l5 " w 15 + "ll\ 

2 V w 15 - "l5 + w 15 + ^15 - "15 + "ll "15 ~ "h + "l5 ~ "l5 ~ "l5 ~ "llj ' 

(21) 

We have computed that these two representations are not entangling. Furthermore we 
have computed that all phase functions of both representations are trivial and so all two- 
qubit projective normaliser elements of are linear normaliser elements. This implies that 
the two representations of [120, 5] given above are not projectively entangling. Hence no 
teleportation group in U{2) which has central quotient isomorphic to the dodecahedral 
group is entangling. 



4.5 A dihedral group as base group 

Since the family of dihedral groups comprises an infinite list we approach the study of their 
projectively inequivalent projective representations analytically. 

The Schur multiplier of a group [T3] is key in calculating covering groups. In particular 
when the Schur multiplier of a group is the trivial group then the group is its own covering 
group. The Schur multiplier M(í?2n) of the dihedral group Dm is given in [Ï9] as 

M(D 2n ) = Z gcd{2)n) . (22) 

This splits the analysis of the dihedral group Z?2n hito the case of odd n where D<m covers 
itself and even n where we obtain the so called binary dihedral groups as covering groups. 
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4.5.1 when n is odd 



We claim that there are no entangling teleportation groups with central quotient isomorphic 
to Z?2n when n is odd. 

As noted above, D2 n , for odd n, is its own covering group. The only representations we 
need to consider in looking for entangling teleportation groups are the the irreducible linear 
ones of Z?2n- When Dm is presented as 

D 2n = (a, b\a n = 1, b 2 = 1, bab = a' 1 ) (23) 

we find, from [TH], the r th irreducible representation p r with < r < ^ may be taken 
to be 

*<•> " (o • * (6) = (! í) ■ <24) 

In particular we are only interested in the faithful representations when gcd(r,n) = 1. 
We prové that the teleportation group G for r = 1 is not projectively entangling and the 
argument for general r follows similarly. Indeed for each r we get the same set of matrices 
(but associated to group elements in different ways) and the normalising property is a 
property of the collection of matrices only £IS cl set. 

Proposition 1. Let G = (a = ^ n J^i ^ ,S = (5q)/ — ^2n for odd n. Then every phase 
function of G takes vàlues in {1,-1}. 

Proof: For any phase function / of G we must have that 

f(A) =<4 and f(B) = (-l) k (25) 

by the orders of A and B respectively. This defines the value of / on all elements of G. It 
is easily verified that if j = and k = 1 then / defines a vàlid phase function of G but that 
for j ^ the group {f(A)A, f(B)B) is never isomorphic to G. □ 

Since every phase function takes vàlues in {1,-1} we must now show that G' is not 
entangling where G' is generated by adding the element ( 7) 1 -i ) to the generators of G . 
So we have 

Cf-(c-(lr .-d.W? (26) 



, u£J ' \1 0, , 
Proposition 2. G 1 is not entangling for any odd n > 3. 

The proof of the above is split into three lemmas from which the result follows. 

Lemma 8. Every two-qubit normaliser gate for G' is a generalised permutatíon matrix. 
That is it has exactly one non-zero entry for each row and each column. 

Proof: Let us write a general two-qubit normaliser gate of G' as ./V with matrix entries 
JVjjfe. We write C\ = C (g> I and C 2 = I <8> C with C as in eq. f[26|) . For any j G {0, 1, 2, 3} 
we have 

{NdN% = ^(lA^-il 2 + |iV i2 | 2 ) + u£(\N j3 \ 2 + |iV,4| 2 ). (27) 
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Since NCiN^ € G' each entry (NC\N^)j must be zero or a power of u> 2n and since n > 3 
we must have either 

lA^-il 2 + lA^I 2 = or \N j3 \ 2 + \N j4 \ 2 = 0. (28) 

Similarly 

(NC 2 N% = uj 2n (\N n \ 2 + |7V i3 | 2 ) + v£(\N j2 \ 2 + \N j4 \ 2 ) (29) 
from which we conclude that either 

l^'il 2 + l^isl 2 = or |iV i2 | 2 + \N j4 \ 2 = 0. (30) 

From equations [28] and [30] we see that N must have exactly one non-zero entry per row and 
since N must be unitary we conclude it is a generalised permutation matrix.D 

Lemma 9. If there exists an entanglíng two-qubit normaliser gate for G' then there also 
exists a diagonal entanglíng two-qubit normaliser gate for G' . 

Proof: Suppose that N is any entangling 2-qubit normaliser for G' . From lemma [8] we 
have N = DP where D is diagonal and P is a permutation. The three 2-qubit normaliser 
gates I <S> B, B ® I and B ® B (with B as in eq. f)26fí ) are permutations that interchange 
00 with 01,10 and 11 respectively. Hence if R is a suitably chosen one of these three, then 
we get an entangling 2-qubit normaliser N' = NR = DP' where P' is a permutation that 
leaves 00 fixed. If S denotes the swap gate (which is a 2-qubit normaliser for any group) 
and Cx denotes the controlled NOT gate, then the six possible choices of P' can be written 
as I, S, C X ,SC X , C X S and SC X S. If P' = I then N' is diagonal. If P' = S then SN' is 
a diagonal entangling normaliser. For P' = Cx recali C 2 = I <g> C (with C € G' as in eq. 
(|26p ). Then a direct calculation shows 

N'C 2 N'ï = diag(Lü 2n , , , cu 2n ) i G' . (31) 

Hence we cannot have P' = Cx- Similarly if P' were SCx,CxS or SCxS we could pre- 
and/or post-multiply N' by S to obtain a normaliser again of the form N" = D'Cx with 
D' diagonal. Hence these three cases of P' are also excluded and in all allowed cases, the 
existence of N implies the existence of a diagonal 2-qubit entangling normaliser. □ 

Lemma 10. No diagonal two-qubit normaliser gates for G' are entanglíng. 

Proof: Let us take an arbitrary diagonal two-qubit normaliser gate of G' which we may 
write up to phase as D = diag(l,a,b,c). We then see that for DB\D^ to be in G' there 
must be j, k such that 

DB ^ = B ) = ® ^ = I . .-ï» ü n I • (32) 

\ ca 

Eliminating b from the above we see that since n is odd we must have k = or k = n. 
Similarly from DB 2 D^ we have l, m such that 
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and eliminating a gives l = or l = n. By solving for a, b, c in the four cases of k, l = 0, n 
we see that in each case D can be written as a tensor product of two gates and hence is not 
entangling. □ 

We have now completed the proof of proposition [2j This completes the result that no 
dihedral group Z?2n when n is odd forms the base group of an entangling teleportation 
group. 

4.5.2 when n is even 

We claim that every dihedral group D2 n , where n is even, is isomorphic to the central 
quotient of an entangling teleportation group. Indeed introduce m defined by n = 2m and 

Z l/m = (l ° ) G m = {X,Z x l m ). 

\ U w 2m / 

A straightforward calculation shows that p defined by p(a) = Z x l m and p(b) = X provides 
a projective representation of Z>2n as presented in eq. ([23]) and G m /Z(G m ) = P>i n . 

For m = 1 the projective normalisers are clearly just those of the Pauli group given in 
lemma [TJ 

For m > 2 we find that Z l / 2m = diag(l,u;4 m ) is a normaliser of G m . This is the 
generalisation of the phase gate P from the Pauli group normaliser (and now H is no longer 
a normaliser). Also CZ is a normaliser of G m ® G m so G m is entangling for all m G N. Our 
computer àlgebra procedures for small m vàlues showed that all normalisers of G m (g> G m 
are generated from CZ with SWAP and G m (8) G m included. 

Finally we show that there are no other teleportation groups G C U(2) with central quo- 
tient Z?2n (n even) that are not unitary equivalents of central extensions of G m above. Indeed 
the covering group of Z?2n (n even) is the binary dihedral group Q^n with presentation[20j 

Q 4n = { a ,b: a 2n = 1, b 2 = a n , b^ab = a' 1 } 

and a complete set of faithful irreducible representations on C 2 is given bv|20j 

Pr(a)=M" ^r) , Pr(6) = (J with gcd(r, 2n) = 1. 

The case r = 1 reproduces G m above and for each r, p r comprises the same set of matrices 
(up to phase múltiples) so we get no new projective normaliser structures as r varies. This 
completes the proof of theorem [2j 

5 Conclusions 

We have identified all entangling teleportation groups in U(2) and seen that they comprise 
only a mild generalisation of the Standard qubit Pauli group. Also the associated projective 
normalisers, apart from extra roots of Z, are already present for the Pauli group case. Thus 
the qubit case appears to be of limited scope in generating new classes of classically simu- 
latable circuits, but there are yet further cases and gener alisations worthy of investigation 
which we list as open qüestions. 



18 



Firstly we may consider higher vàlues of d. The Pauli groups may be naturally gen- 
eralised to arbitrary dimension d and the associated normaliser groups were analytically 
characterised for all prime d in [2T]. In this case it is found that an analogue of lemma 
[T] holds (with H, P and CZ being replaced by suitable one and two qudit operations as 
given in |_21j). We applied our computational programs to a further few chosen examples 
of teleportation groups in Í7(3) but did not find any further interesting entangling ones. 
We were unable to exhaustively treat all base groups of Í7(3) because of the increased size 
of the groups involved. Thus it would be advantageous to further develop the study of 
computational procedures for projective normalisers, inventing algorithms that search over 
more restricted spaces of vàlues. 

We have considered the projective normaliser structure of G and G®G only for matrix 
groups G that act irreducibly. But more generally if a matrix group acts reducibly it is not 
clear how its normaliser structure relates to that of its irreducible parts. This may provide 
an avenue for constructing further interesting examples of entangling normalisers for groups 
acting in dimensions d > 3. (For the case d = 2 that we have considered exhaustively any 
reducible group is diagonal). 

Our exhaustive qubit results indicate that the Pauli matrix group is very special in 
possessing a suitably rich variety of normalisers. In this vein it would be particularly 
interesting to identify a mathematical "signature" property of a given matrix group G 
whose validity signals the existence of non-trivial projective normalisers. 

Finally we point out that lemma [T]asserts a remarkable structural property of normalis- 
ers of the Pauli group V and its tensor powers J>® n viz. that for levels n > 3 there are no 
new normalisers beyond those generated from circuits of n = 1 and 2 normalisers. (This is 
also true of the generalised Pauli groups in prime dimension [21]). The proof of this property 
utilises many extra properties special to the Pauli matrices. Thus we may ask: are there 
teleportation groups G (even in dimension d = 2) such that G® n for n > 3 has normalis- 
ers that are not expressible as composites of n = 1,2 normalisers? In our computational 
analyses we have considered projective normalisers only up to tensor square G <8> G and 
it remains open whether or not there may exist entangling normalisers for n > 3, even in 
the case that they are absent for n = 2. Any mathematical signature property of the kind 
mentioned above would be helpful in addressing this fundamental issue. 
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